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Abstract. Coupled aggregation and sedimentation processes were studied by means of three dimensional 
computer simulations. For this purpose, a large prism with no periodic boundary conditions for the sedi- 
mentation direction was considered. Furthermore, three equally sized and mutually excluded regions were 
defined inside the prism, a top, a middle and a bottom region. This allows to study the time evolution 
of the cluster size distribution and the cluster structure separately for each region. The mass distribution 
profile and the center of mass position were also accessed as a function of time. For the bottom region, 
the effects of the sediment formation on the kinetics of growth and on the cluster structure were clearly 
observed. The obtained results not only agree with the experimental data obtained by AUain et. al. |^] and 
with the simulations made by Gonzalez M but also allow to gain further insight in the details. 
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1 Introduction 

To a larger or smaller extent, all real aggregation processes 
are influenced by the presence of an external gravitational 
' field. Hence, the aggregation phenomena are usually fol- 
, lowed by sedimentation or, in many cases, both are present 

■ simultaneously. Some natural examples are the delta for- 
mation at a river estuary and the settling of bacteria clus- 

. ters in quiet water jsj . As technological examples one may 
' cite water treatment for human consumption (clarifying) , 
, effluents treatment, and a large number of precipitation 

■ techniques employed by the chemical industry 

[ Pure irreversible aggregation processes have been well 
I studied and are described in the literature. The first equa- 
tion for describing the aggregation kinetics of diluted sys- 
tems was given by Smoluchowski in the early 1900s 1^,^. 
This equation defines an infinite two-dimensional matrix 
of kinetic rate constants, known as kernel, which accounts 
for the physical characteristics of the system. The time 
evolution of the cluster size distribution is obtained by 
solving the Smoluchowski's rate equation for a given ker- 
nel. On the other hand, much information about the aggre- 
gation processes is also contained in the clusters structure. 
Here, Smoluchowski's rate equation does not provide use- 
ful information and hence, other techniques such as light 
scattering experiments and simulations are required for its 



study. Significant contributions in these fields are provided 
by Lin, Weitz and coworkers (light scattering experiments) 



and Meakin, Family and Gonzalez (simulations) 
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Pure sedimentation phenomena have also been well 
studied 
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Here, an important issue is the variation 
of the settling velocity with the volume fraction |p^ , p!7| . 
When the concentrations become large enough, the re- 
verse fiow of fiuid necessary to compensate the volumet- 
ric flow of particles plus the associated fiuid contribute 
to decrease the Stokes velocity (this is known as back- 
fiow effect). Furthermore, since the volume fraction be- 
comes larger when moving to the fiask bottom, an os- 
motic pressure appears opposed to the gravitational field. 
Finally, for larger volume fractions the inter-particle dis- 
tances shorten and hence, the role of hydrodynamic forces 
and inter-particle interactions become important. 

Although aggregation and sedimentation phenomena 
are closely connected, there is not an extended literature 
dealing with them simultaneously. This is, at least par- 
tially, due to the mathematical difficulties, which appear 
when following a formal analytical treatment |l^. Con- 
sequently, computer simulations become a useful alter- 
native tool for studying and predicting the behavior of 
real aggregation-sedimentation systems. Recently, Gonza- 
lez and Leyvraz worked on computer simulations for elu- 
cidating the experimental results found by AUain et. al. 
[0,|l|]. They found that increasing the intensity of the ex- 
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ternal field leads to an increase in the cluster fractal di- 
mension. Their model, however, did not consider cluster 
deposition (only a sedimentation velocity was added to 
the Brownian motion) [^ . Later, the simulations were 
improved by considering a cubic lattice of size L, with pe- 
riodic boundary conditions in the three spatial directions 
(as in the previous model), but now taking the clusters 
out of the cubic box with a probability related to the to- 
tal distance that they had moved downwards (see ref. 
for details). Although this model does not consider the 
mass distribution dependence on height, it was useful for 
explaining an increase of the initial particle concentration 
required to have gelation when sedimentation is present. 

In this paper, we further study the coupled aggrega- 
tion and sedimentation phenomena by simulations, but 
now considering a large prism with no periodical boundary 
conditions for the sedimentation direction (this requires a 
macroscopic prism height). For the other two horizontal 
directions, we still impose periodical boundary conditions 
in order to represent, in an average way, a portion of the 
whole system. This allows us to study not only average 
quantities but also their change as a function of height. 
Furthermore, it was possible to observe and study the sed- 
iment formation. 



2 Simulation method 

The simulation processes were carried out off lattice on a 
square section prism of side L and height H. Inside, A^o 
identical hard spheres of radius a were randomly placed 
avoiding overlapping among them. Since aggregation and 
sedimentation processes are simulated, there are two con- 
tributions to the movement. One corresponds to the Brow- 
nian motion and the other to the Stokes sedimentation 
velocity. As the simulations were done for very dilute sys- 
tems, backflow and hydrodynamic forces were not taken 
into account. This assumption will have to be checked 
since local high concentrations may appear somewhere in 
the prism. The time step was fixed by the relationship 



(1) 



where Di = ksT / {Qirrjo) is the monomer diffusion coeffi- 
cient, ksT is the thermal energy, t] is the solvent viscosity 
and /s is the Brownian length step. Hence, monomers are 
always moved Ib in a random direction plus the Stokes 
contribution 

Is = vstQ = llPe/Qa (2) 
where vs = 2{p — p^jga? / (^rf) is the monomeric Stokes 
velocity, p is the particle density, po is the fluid density, 
g is the earth gravitational constant and Pe — Ana"^ [p — 
Po)g/kBT is the Peclet number. For i-sized aggregates, 
we assume the relationship Vg = ai^/'^^ for the radius 
of gyration and therefore the diffusion coefficient Di = 
ksTi'^/''-^ /(Gnrja) and the Stokes velocity vs = 2(p — 



po)ga^i 



' /9r]. Here, it was implicitly assumed that the 



hydrodynamic fractal dimension and the cluster fractal di- 
mension are equal. Now, in order to account for the Brow- 
nian motion, we move the aggregates Ib in a random direc- 
tion only when a random number, ^, uniformly distributed 



in [0,1], is less than the ratio between the aggregate diffu- 
sion coefficient and the monomer diffusion coefficient, i. e. 
when ^ < . In case that the Brownian movement is 

refused, the corresponding Stokes contribution 



Is = vsto 



/6a 



(3) 



is accumulated in a memory place, associated with the 
particular cluster, for being considered in the following 
time intervals. Once a Brownian contribution is accepted, 
the accumulated Stokes contribution is added to the clus- 
ter motion and the corresponding accumulated memory is 
reset to zero. Furthermore, if a given accumulated Stokes 
contribution exceeds a Brownian step, Ib, then the ag- 
gregate is moved downwards and the accumulated mem- 
ory is also reset to zero, thus avoiding a large total step 
length. The algorithm minimizes the times that clusters 
are moved, which is the most time consuming contribu- 
tion to the execution time. After any movement, the re- 
gional configuration is checked for overlaps. In case that 
any overlap is found, it is corrected by placing the recently 
moved cluster in touch with the other to form a new ag- 
gregate, i. e. diffusion limited cluster aggregation (DLCA) 
conditions are imposed. 

Periodic boundary conditions were established for the 
two non sedimentation directions, x and y. Hence, the 
system may be understood as a small portion of a macro- 
scopic one. However, for the sedimentation direction, z 
(we define the prism bottom as z = and the prism top 
as z = 1), no periodic boundary condition was imposed in 
order to naturally obtain a change in properties with the 
prism height and a cluster deposit at the prism bottom. 
This forces us to use a high prism since the height of the 
macroscopic system coincides with its vertical length. This 
could be done while keeping a reasonable total number of 
particles, i. e. working with diluted systems. As the sys- 
tem symmetry was broken for the z direction, we look for 
this effect on the cluster populations and on the cluster 
structures. 

For this purpose, the system volume was arbitrarily 
divided in three equally sized and mutually excluded re- 
gions: a top, a middle and a bottom region. Hence, we 
may define the weight average cluster size for each region 
as 



(4) 



where the symbol \h — top, middle or bottom refers to 
the corresponding region and Ui is the number of i-sized 
clusters. Furthermore, the vertical position of the system 
center of mass is also defined as 



(5) 



where i; is the size of cluster I and zf™ is the vertical 
position of the center of mass of cluster I. 
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Fig. 1. Three dimensional representation of a simulated sys- 
tem settling under Pe = 0.1. The prism dimensions are H — 
1500a and L = 250a (H = 1500a was established instead of 
H — IQQQOa for an easier representation). The images were 
captured at times 0.1 s a), 112 s b), 630 s c) and 1120 s d). 



The cluster structure is studied by calculating the fol- 
lowing radii of gyration every time a new cluster is formed 



gx 



(6) 



^E^(i-i.ri) 



where is the distance between the particle i and the 
cluster center of mass and x^, yi and Zi are its compo- 
nents. Those radii of gyration allow us to evaluate the 
average ratios (rgx) / {rgz)\h and {rgy)/{rg^)\h and to ob- 
tain three cluster fractal dimensions from rg\h = ai^^'^f^'^ 
(considering only clusters containing more than 15 parti- 
cles) . The former quantities account for the average shape 
of the clusters, i. e. {rgx)/{rgz) > 1 indicates that the 
clusters tend to be elongated in the z direction and, on 
the contrary, (rgx) / (rgz) < 1 points to wider structures. 
{rgx ) I {vgz ) and {vgy ) I {vgz ) shouM be similar but not equal 
due to statistical fluctuations and so, their difference rep- 
resents an estimation of their uncertainty. 

It should be noted that the cluster motion depends on 
the fractal dimension, which is not known a priori. This 
forces us to estimate df and to iterate the simulations 
whenever the introduced df differs from the obtained df 
in more than 0.05. 

The parameters employed for the simulations arc the 
following: a monomer radius a = 315 nrn, a Brownian 
step length Ib = a/ 2, a prism height H = 10000a, a sec- 
tion side L = 250a and a particle volume fraction of (/> = 
6.7 X 10~^ (iVo = 10000). The particles were considered to 
be dispersed in water at 20° C. Different runs were per- 
formed for different Peclet numbers. In real experiments, 
a large change of the Peclet number may be achieved by 



changing the strength of the external field by centrifuga- 
tion. Another possibility for realizing small changes of the 
Peclet number consists in using particles with different 
densities. 



3 Results 

3.1 Overview 

Figure ^ was constructed to give an overview of the cou- 
pled aggregation and sedimentation processes. It shows a 
three dimensional representation of a system settling un- 
der Pe = 0.1. The images were captured from simulations 
at times 0.1 s a), 112 s b), 630 s c) and 1120 s d). In order 
to easily represent the images, here we set H — 1500a in- 
stead oi H = 10000a, which is the value employed for the 
study. As can be seen for the beginning of the process, 
figure |l] a) shows a uniformly distributed and monodis- 
perse system. After a few seconds, figure |^ b), the sys- 
tem aggregates forming a wide cluster size distribution. As 
imposed, the larger clusters move faster downwards and 
hence, they are mostly seen close to the prism bottom. As 
the times goes on, those clusters start forming the sedi- 
ment while the dispersion clarifies by loosing mass. The 
first settling aggregates that arrive to the prism bottom 
continue their Brownian motion although they are not al- 
lowed to move further downwards. Hence, their movement 
is almost restricted to two dimensions, since their weight 
makes difficult their upward motion. Consequently, they 
collide among each other forming the sediment. In the fol- 
lowing time intervals, the system looks like the one shown 
in figure |l| c). Here, it is observed that the mass concen- 
tration depends on the distance from the bottom of the 
prism. At the top, the system is almost clear while the 
dispersion looks more concentrated at the bottom. Fur- 
thermore, a large sediment is also observed. It gains mass 
from every closed enough settling cluster. On the other 
hand, the number of clusters remaining in the dispersion 
strongly decreases and the aggregation rate of the dis- 
persion slows down. Consequently, aggregation becomes 
highly improbable for those small remaining clusters and 
a long time is needed for their deposition. This situation 
is clearly seen in figure l| d) . 

It should be pointed out that this first overview is in 
good agreement with the experimental observations made 
by AUain et. al. for calcium carbonate colloidal suspen- 
sions |jl[. Their observations were performed experimen- 
tally in a 800 mm high cell. Furthermore, they emphasize 
that the different phases of settling cannot be identified 
separately for smaller cells. This is not surprising since 
their experiments were characterized by a Peclet number 
close to lO"'^. 



3.2 Further inside in the details 

Figure ^ a) shows the time evolution of the cluster-size 
distribution for the upper portion of the system, h — 
top, aggregating under Pe — 0.1. As expected, the first 
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Fig. 2. a) Time evolution of the normalized cluster-size dis- 
tribution for the upper portion of the system, h = top, and Pe 
= 0.1. The points are grouped in logarithmically spaced inter- 
vals ((□) monomers, (o) 2 and 3-mers, (a) 4 to 8-mers, (v) 9 
to 18-mers, (o) 19 to 38-mers and (+) 39 to 88-mers. b) The 
corresponding weight average cluster size, n™. 
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Fig. 3. a) Time evolution of the normalized cluster-size distri- 
bution for the middle portion of the system, h = middle and 
Pe = 0.1. The points are grouped in logarithmically spaced 
intervals ((□) monomers, (o) 2 and 3-mers, (a) 4 to 8-mers, 
(v) 9 to 18-mers, (o) 19 to 38-mers, (+) 39 to 88-mers and (x) 
89 to 200-mers). b) The corresponding weight average cluster 



stage of the aggregation process evolves in time similarly 
to the DLCA regime, i. e. the sedimentation effects are not 
yet important. Once the oligomer concentration becomes 
higher, the subsystem starts to loose mass due to sedimen- 
tation and so, the formation of larger clusters is biased. 
Furthermore, the characteristic bell shape evolutions of 
oligomers narrow since they disappear not only by reac- 
tion but also by leaving the subsystem. When comparing 
the obtained time evolution of the monomer concentration 
with the one corresponding to the DLCA regime, differ- 
ences are also observed for the later stages. For the pure 
DLCA regime, the monomer concentration decays faster 
than in the case where sedimentation effects are consid- 
ered. The reason for this is that a larger cluster concen- 
tration remains in the subsystem and hence, a higher ag- 
gregation velocity is achieved. Moreover, since for DLCA 
processes the monomeric aggregation rate constants are 
largest, their curve crosses the curves for the larger clus- 



ters leading to a smaller monomer concentration P5[ . This 
is not observed when sedimentation is present. 

The corresponding weight average cluster for 
h — top and Pe = 0.1 is shown in figure || b). Again, at 
the beginning of the processes n^, evolves in time similarly 
to the DLCA regime. As time goes on, two effects tends 
to change its evolution. On the one hand, the different 
sedimentation velocities of clusters leads to an increase in 
their collision frequency and consequently, the average riw 
tends to increase its rate of change. On the other hand, 
the larger aggregates exit the subsystem more often than 
the smaller ones tending to decrease the average n^. For 
a given time, the latter effect prevails over the former and 
so, the average cluster size peaks. For larger times, this 
situation makes the average to monotonically decrease 
until monomers are the only species in the subsystem. 
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Fig. 4. a) Time evolution of the normalized cluster-size distri- 
bution for the bottom portion of the system, h — bottom and 
Pe = 0.1. The points are grouped in logarithmically spaced 
intervals ((□) monomers, (o) 2 and 3-mers, (a) 4 to 8-mers, 
(v) 9 to 18-mers, (o) 19 to 38-mers, (+) 39 to 88-mers and (x) 
89 to 200-mers). b) The corresponding weight average cluster 
size, n„. 



The time evolution of the cluster size distribution and 
the corresponding average, n^, for the second region, h 
= middle, are shown in figure ra. Although the evolutions 
are similar than those shown for the upper region of the 
prism, some differences are found. Since this subsystem 
also gain mass from the upper region while loosing it to- 
wards the bottom, it is capable to develop larger struc- 
tures and to hold a large number of clusters. This is clearly 
seen in figure ^ a) where a large number of 39 to 88-mers is 
shown and even a few 88 to 200-mers appear. Furthermore, 
the weight average cluster size peaks for larger times and 
reaches a higher value (compare figure || b) with |^ b)). Fi- 
nally, it is observed that the monomer population becomes 
almost constant at approximately t — 4000 s, which is the 
same time at which the upper subsystem has almost lost 
the remaining monomers. This fact is explained by means 
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Fig. 5. Mass distribution profile, m(z), along the sedimenta- 
tion direction, z, represented in log-normal axes for different 
times. The data were obtained for Pe = 0.1. The symbols □, 
o, A and V represent the mass distribution for t = 112 s, 355 
s, 1122 s and 3550 s, respectively, z = corresponds to the 
prism bottom and 2 = 1 to the top. The dashed lines indicates 
the boundaries of the three defined regions. 



of the monomeric transfer from the upper region towards 
the middle subsystem. 

Figure ^ shows the time evolution of the cluster size 
distribution and weight average cluster size for the bottom 
region of the prism. Since this subsystem gains mass from 
the upper regions, its particle concentration continuously 
increases and hence, some remarkable differences appear. 
The most important is that the average n„ becomes a 
monotonously increasing function of time. Moreover, nu, 
seems to follow an asymptotical power behavior, n^, ~ t'^, 
where fc ~ 4. Due to the finite extension of the studied 
system, the asymptotical behavior cannot be prolonged in 
time and hence, n^, diminishes its rate of change until the 
remaining clusters settle and become part of the sediment. 
Another difference is that larger clusters are formed at the 
bottom faster than in the bulk. Furthermore, they also 
react faster to form the sediment and so, their time evolu- 
tion shows a double peak. The former peak (the smallest) 
corresponds to those clusters aggregating at the bottom 
to form the sediment and the second to the bulk clus- 
ters. Monomers also behave differently. During the final 
stages their concentration abruptly change from diminish- 
ing to increasing. This is a consequence of the monomer 
flow coming from the middle region that prevails over the 
small remaining aggregation rate. For even longer times, 
monomers eventually collide with the sediment and hence, 
they finally disappear. 

In order to study the mass distribution profile along 
the sedimentation direction, z, the system was subdivided 
in 20 slices. For each one and for a given time, its total 
mass was calculated by m(z) = in\^ and represented 
as a point in figure |. This fi gure was constructed for Pe 
~ 0.1, i. e. corresponds to the data shown in figures |^, 
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Fig. 6. Time evolution of the weight average cluster size, n^, 
for different Peclet numbers, Pe, and for the upper a), middle 
b), and bottom c), portions of the system. The symbols □, o, 
A, V and o correspond to the evolutions obtained under Pe = 
1, 0.1, 0.01, 0.001 and 0.0001, respectively. The dashed lines 
represent the asymptotical DLCA behavior. In figure c), the 
continuous lines are drawn as a guide to the eye. 



^ and 1^. It is observed for time 112 s that the mass dis- 
tribution is practically uniform. Only for the upper and 
bottom slices, slight deviations from the mean value are 
observed. While the system evolves in time, those devia- 
tions becomes larger due to the settling process. At 355 s, 
the mass concentration of the upper slice is about one half 
of the initial value and the mass concentration of the bot- 
tom slice is approx. the double. The other slices, however, 
keep their initial mass concentration. After a sufficiently 
long time (see 1122 s), no remaining particles are observed 
in the upper slice and the mass concentration diminishes 
strongly in the two contiguous slices. Nevertheless, the 
other slices do not change their mass concentration ex- 
cept for the bottom one, which contains the sediment. 
This situation changes for very long times, where the ini- 
tial mass concentration is not maintained any more in any 
slice of the system. This can be seen for time 3550 s. Here, 
the sediment contains more than half of the total mass. It 
should be noted that the mass concentration is not larger 
than the initial concentration except for the bottom slice. 
This means that the assumption made in section |^ of ne- 
glecting the backflow effect and the hydrodynamic forces 
due to the low concentration is valid for almost the whole 
system. 

The time evolution of the weight average cluster size, 
n^,, for the Peclet numbers Pe — 1, 0.1, 0.01, 0.001 and 
0.0001 is shown in figure ^. Again, three different plots 
are employed for representing the evolution of each re- 
gion. In addition, all figures include a unity sloped straight 
line, which represents the asymptotic evolution of Uw for 
DLCA. As expected, the sedimentation effects appears at 
shorter times as the Peclet number is increased. This is 
clearly seen for the three prism regions. On the other hand, 
the beginning of the coupled aggregation and sedimen- 
tation process evolves similarly to one following a pure 
DLCA regime for all the studied Peclet numbers. This in- 
dicates that a very large Pe is needed for appreciating the 
sedimentation effects at the beginning of the processes. 
Figures ^ a) and b) show that the DLCA asymptotic evo- 
lution is, generally, not surpassed by the curves that con- 
sider sedimentation. In fact, when settling effects appear 
the weight average cluster size starts diminishing. Further- 
more, it is observed for all Peclet numbers that the 
curves peak earlier and reach lower values for the upper 
region. On the contrary, in figure g c) it is observed that 
the curves behaves completely diuerent. Here, the sedi- 
mentation effects also change the evolution by strongly 
increasing its slope. We found for Pe = 1, 0.1 and 0.01 a 
practically constant slope of approximately 4, which is, by 
far, larger than the DLCA asymptotic slope. The continu- 
ous lines shown in figure ^ c) for this Peclet numbers were 
drawn by considering this slope. For Pe = 0.001 a slope of 
approx. 2.5 was obtained and for Pe — 0.0001 no signif- 
icant deviation from the pure DLCA slope was observed. 
It is likely that, at least for our simulation conditions, a 
limiting value of 4 is achieved for the asymptotic slope. 

The vertical position of the system center of mass, Zcm, 
is plotted against the dimensionless time 9 = t{Pe)°' /tag 
in figure 0, for the set of Peclet numbers. Here, tag = 
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Fig. 7. a) Zcm as a function of S for Pe = 1 (□), 0.1 (o), 0.01 
(a), 0.001 (v) and 0.0001 (o), represented in a semilog plot. 



2/Coi^fi™°' = 352s is the characteristic aggregation time, 



Co = 5.12 X 10^* m~'^ is the initial particle concentra- 
tion, iff™°' = SksT/ (3ry) is the dimer formation Smolu- 
chowski's rate constant, which takes the value 11.1 x 
lO^^^rn^s^^ for the assumed solvent conditions. In ad- 
dition, the parameter a was introduced in order to at- 
tempt to define a single master curve for the complete 
set of curves. As figure ^ clearly shows, the master curve 
is yielded for a = 0.74. This curve takes a value closed 
to 0.5 at the beginning of the processes since the mass 
is randomly scattered at time zero for each system, Zcm- 
At first stages and for all Pe, Zcm does practically not 
evolve in time. Once the aggregation process leads to large 
enough clusters, the settling process starts taking place. 
This starting point strongly depends on the Peclet num- 
ber. The larger the Peclet number is, the sooner the set- 
tling process appears since a smaller cluster size is needed 
for obtaining an appreciably sedimentation velocity. Once 
the settling effects appear, the system center of mass moves 
quickly towards the prism base until the mass transfer 
from dispersion to sediment is finished. 

As explained in section ^, the cluster fractal dimension 
was assessed by means of the radii of gyration method and 
calculated for each defined region. For Pe = 1 and for the 
upper regions, it was not possible to obtain a reliable frac- 
tal dimension value since not enough large clusters were 
formed in these regions. For obtaining a reliable statis- 
tics, three simulation runs were performed for each Peclet 
number. Pure DLCA simulations were also carried out as 
a reference. The results are shown in table |l|. For Pe — 0, 
i. e. pure DLCA, the values are consistent with the ones 
reported in the literature |21,2^. Furthermore, no signifi- 
cant difference appears among the regions, as expected for 
a regime where no sedimentation velocity was imposed. 
For larger values of the Peclet number it is observed that 
the cluster fractal dimension increase for the top and mid- 
dle regions. This is in agreement with the simulations re- 
sults obtained by Gonzalez, in which no extra effects are 
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Fig. 8. The cluster size, i, as a function of the normalized 
radius of gyration, rg/a, for Pe = 0.01 and for the middle 
region. The solid lines represents the linear fits obtained for 
the ranges rg/a e [2,6] and rg/a e [6, 18], respectively. 



needed, such us cluster restructuring, for explaining the 
increase of the fractal dimension Nevertheless, Gon- 
zalez obtained df = 2.27 for the larger structures and for 
Pe — 0.01, which is larger than our results. This may be a 
consequence of the higher concentrations employed in his 
simulations as discussed further in the text. In addition, 
the dflmiddie valucs seems to be slightly larger than the 
df\top values. This may be due to the average downward 
distance that the corresponding clusters moved. Since this 
distance is longer for the clusters inside the middle region, 
the effect of increasing the fractal dimension due to the 
Stokes velocity becomes more evident. 

For Pe = 0.01 and for the middle region, figure |8| shows 
the cluster size, as a function of the normalized radius of 
gyration, rg/a. The data were fitted over two ranges and 
different slopes were obtained. This means that a crossover 
between the two different power lows may, at least, not be 
discarded. This behavior was already reported in ref. 
where it is explained that smaller clusters show a roughly 
DLCA behavior, whereas larger ones behave differently 
since they are more influenced by the external field. Fur- 
thermore, it is expected that the crossover cluster size de- 
creases as Pe increases as reported in ref. This is not 
easy to see in our simulations since the clusters disappear 
from the middle region before growing enough, and this 
leads to poorer statistic results. It should be pointed out 
that table |l| only reports average values (considering only 
those clusters containing more than 15 particles), which 
were employed to dictate the cluster motion, i. e. no dis- 
tinction was made between the fractal dimensions of small 
and large clusters. 

The fractal dimensions obtained for the bottom region 
varies in a quite different way with the Peclet number than 
for the upper regions. Here, it is observed that df strongly 
decreases as the Peclet number increases. It was expected, 
however, to obtain even larger fractal dimensions than the 
ones obtained in the upper regions. The fact of obtain- 
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Table 1. Cluster fractal dimensions, and average ratios {rgxy)/{rgz)\h as a function of the Peclet number, h refers to the 

top, the middle and the bottom regions. The last column indicates whether the system percolates at the prism base. 



Pe 


df\top 


df 1 middle 


df\bottom 


{'''g^y) / {''ffz) I'op 


{Tgxy)/ {rgz)\middle 


(j'gxy) / {rgz)\bottom 


percolation 


1 






1.67 ±0.05 


1.01 ±0.01 


0.99 ±0.03 


0.90 ± 0.02 


Yes 


0.1 


1.88 ±0.07 


1.89 ±0.06 


1.74 ±0.05 


1.00 ±0.05 


1.01 ±0.02 


0.95 ±0.02 


Yes 


0.01 


1.80 ±0.06 


1.83 ±0.06 


1.82 ±0.05 


0.99 ±0.05 


1.00 ±0.02 


0.99 ± 0.02 


Yes 


0.001 


1.82 ±0.06 


1.84 ±0.06 


1.81 ±0.05 


1.02 ±0.05 


1.00 ±0.02 


1.01 ±0.02 


No 


0.0001 


1.83 ±0.06 


1.80 ±0.06 


1.82 ±0.05 


1.00 ± 0.05 


0.99 ±0.02 


1.00 ±0.02 


No 





1.78 ±0.06 


1.74 ±0.06 


1.76 ±0.05 


0.99 ± 0.05 


1.00 ±0.02 


1.00 ±0.02 


No 



ing a fractal dimension as low as df — 1.68 ± 0.05 for 
Pe = 1 clearly indicates that the sediment formation is 
changing in some way the growth process. Since the clus- 
ter movement is restricted almost to two dimensions at 
the prism base and since the fractal dimension yield for 
DLCA processes at two dimensions is df — 1.45 ± 0.05, 
it is not very surprising to find cluster structures char- 
acterized by fractal dimensions ranging between 1.75 and 
1.45 1^,^. This indicates that for increasing values of 
the Peclet number, the average cluster size of the settling 
aggregates that arrive to the prism bottom decreases and 
hence, most clustering reactions take place in the bottom 
region. 

The average ratio {rgxy) / {rgz)\h = {{rgx) / {rgz)\h + 
{'''gy) /(''ffz)U)/2 obtained for different Peclet numbers is 
shown in table |[ As explained in section ||, this quan- 
tity accounts for the average shape of the clusters. In case 
that {rgxy) / {vgz) > 1, the clusters tend to be elongated 
in the z direction and when (rgxy) / (rgz) < 1 the clus- 
ters are shorter. For the upper regions, no significative 
deviation from unity is observed not even for the higher 
Peclet numbers. This indicates that the clusters do not 
show a preferential growth direction. This is in good agree- 
ment with AUain et. al. experimental observations U|l^. 
Nevertheless, our results do not agree with Gonzalez find- 
ings For Peclet numbers as high as 0.1 and for much 
more concentrated systems (0 = 0.01), Gonzalez found 
that clusters grow faster in the sedimentation direction. 
This is very likely to occur since clusters move mostly 
downwards. For diluted systems, on the other hand, when 
the Peclet number is as high as 0.1, the clusters settle so 
fast that they have not the opportunity to grow before 
arriving to the prism base. This fact explains why we do 
not obtain {vgxy) / {rgz) significatively larger than unity 
for the upper regions. The same argument is also valid 
to explain the differences between the fractal dimensions 
obtained by Gonzalez and by us. For the lower region and 
where the sediment grows, for the highest Peclet numbers 
we obtain (rgxy) / (rgz) lower than unity. This is a conse- 
quence of the fact that aggregation is taking place mostly 
at the prism bottom, at which the z direction motion is 
restricted. Hence, the sediment grows covering the base 
of the prism producing local two dimensional percolation. 
This kind of percolation, as shown in table |^, occurs only 
for the largest values of the Peclet number. 

Finally, figure ^ shows the sediments obtained for dif- 
ferent Peclet numbers. It can clearly be seen that for in- 
creasing Pe the sediments become more extended in the 



base plane of the prism and show a higher degree of com- 
pactness. For instance, sediments a) and b) almost cover 
the prism bottom completely. In order to study the struc- 
ture of the sediment, a fractal analysis was realized for 
each sediment cluster. For this purpose, the distances be- 
tween the centers of each particle and the center of a given 
reference particle are calculated. This allows to build a 
function i(r), where i denotes the number of particles 
inside a sphere defined by the reference particle center 
and the radius r. The procedure is repeated by changing 
the reference particle until all particles had been consid- 
ered as reference, while averaging all obtained i{r) func- 
tions. Then, df may by obtained from the relationship 
r ^ i^/'^f . For improving the statistics, three sediments 
were analyzed for each Pe and the corresponding df val- 
ues were averaged. This procedure yields well defined df 
for the smaller Peclet numbers, having df — 1.82 ± 0.05, 
1.91 ± 0.06 and 1.90 ± 0.06 for Pe = 0.0001, 0.001 and 
0.01, respectively. This values are all slightly higher than 
the corresponding values shown in table |l], which may be 
due to differences on the applied methods. Anyway, the 
same tendency for increasing Pe is verified. For Pe — 
0.1 and 1, two different power low behaviors were found 
for the i{r) functions. A short range df — 1.99 ± 0.07 
and 1.89 ± 0.07 and a long range df = 1.58 ± 0.07 and 
1.67 ± 0.07 were obtained for Pe = 0.1 and 1, respec- 
tively. This points towards a crossover from an aggrega- 
tion mechanism which takes place mostly in the bulk to 
another aggregation mechanism restricted almost to two 
dimensions, i. e. processes occurring at the prism bottom. 
Figure |l^ shows the averaged i{r/a) function for Pe = 1. 

4 Conclusions 

Coupled aggregation and sedimentation process were sim- 
ulated by considering a large prism with no periodical 
boundary conditions for the sedimentation direction. Three 
equally sized and mutually excluded regions were defined 
for studying their time evolution of the cluster size distri- 
bution. We found that for the upper regions the shape of 
the time evolution of oligomers narrows due to the effect 
of the Stokes velocity. On the contrary, the population 
dependence on time shows a double peak for the bottom 
region as a consequence of the settling clusters that arrive 
at the prism base. The time evolution of the weight aver- 
age cluster size also changes its behavior depending on the 
region. While it shows a peak for the upper regions, it be- 
comes a monotonously increasing function for the bottom 
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Fig. 9. Three-dimensional perspective plot of the sediments 
obtained for Pe = 1 a), 0.1 b), 0.01 c), 0.001 d) and 0.0001 e). 
No sediment was obtained for Pe = 0. 
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Fig. 10. The i{r/a) function averaged for the three sediments 
obtained under Pe = 1. 

region. Furthermore, the limiting value of 4 was obtained 
for its slope for increasing vahics of the Peclet number. 

The cluster structure was also studied by means of 
measuring the cluster radii of gyration. In agrement with 
the experiments performed by AUain et. al., we obtained 
that there is not a preferential growth direction. We also 
found an increase in the fractal dimension for increasing 
Peclet number in the upper regions. For the bottom re- 
gion, however, a decrease in the fractal dimension and 
the preferential growth directions parallel to the prism 
base were found for increasing Peclet number. In addition, 
for the highest Pe, the sediment fractal dimension crosses 
over from high values, corresponding to the smaller radii, 
to values close to 1.45 ± 0.05, for larger radii, which is the 
accepted value for two dimensional DLCA. These facts are 
explained as a consequence of the restricted motion of the 



settling clusters that have already arrived at the prism 
base. Finally, one may conclude that increasing the Peclet 
number causes the system to increase its tendency to pro- 
duce a two dimensional percolation at the prism base. 
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